Theory of bright-state stimulated Raman adiabatic passage 
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We describe analytically and numerically the process of population transfer by stimulated Raman 
adiabatic passage through a bright state when the pulses propagate in a medium. Limitations of the 
adiabaticity are analyzed and interpreted in terms of reshaping of the pulses. We find parameters 
for the pulses for which the population transfer is nearly complete over long distances. 
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I. INTRODUCTION 



A very popular technique for complete population 
transfer in lambda systems uses stimulated Raman pro- 
cesses by adiabatic passage (STIRAP) 0, S H 9- The 
method is based on a so called counterintuitive sequence 
of laser pulses in which the Stokes pulse, coupling the fi- 
nal state and the intermediate excited state, precedes the 
pump, coupling the initial state and the excited state. 
The propagation of a counterintuitive pulse sequence in 
a A system in the STIRAP regime has been theoreti- 
cally investigated in [f|. An efficient population transfer 
has been reported in crystals doped with rare-earth ele- 
ments Pr 3 +:Y 2 Si0 5 0,0, %. We note that, due to their 
large density and scalability, solid media are of signifi- 
cant interest for applications, e.g. in optical data storage 
and processing. Particular solids, e.g. quantum dots, 
colour centers or doped solids combine the advantages 
of atoms in the gas phase (i.e. spectrally narrow tran- 
sitions) and solids (density and scalability). Reference 
[H also shows an alternative efficient method of popu- 
lation transfer where an intuitive sequence of pulses is 
used with a large one-photon detuning. Unlike the STI- 
RAP for which the dynamics adiabatically projects along 
a dark state (i.e. no component of the excited state), 
the intuitive process follows a bright state. It has been 
named bright STIRAP (b-STIRAP). The possibility of 
such a transfer for one atom was predicted and analyzed 
in[i,[T|. 

In this paper we investigate theoretically population 
transfer by b-STIRAP in a medium. A detailed study 
of the system of Maxwell and Schrodinger equations is 
performed. We derive a self-consistent solution of the 
problem taking into account the first order nonadiabatic 
corrections. From the solution obtained we derive the 
conditions of complete population transfer by adiabatic 
passage in a medium at large propagation distances. We 
show that during propagation both pulses experience a 



reshaping. From the obtained solution we derive a crite- 
rion for adiabaticity condition. 

The paper is organized as follows. In Sec. II we ex- 
plain a theoretical model for b-STIRAP. The solutions to 
the propagation equations are presented in Sec. III. In 
Sec. IV we apply the theoretical model to a real experi- 
ment and discuss the results obtained comparing analyt- 
ical and numerical solutions. Section V is devoted to the 
pulse dynamics. The conclusion is presented in Sec. VI. 



II. THE MODEL 

We consider the propagation of near-resonant pump 
and Stokes pulses (of respective frequency to p and u> s ) 
in a medium of A-type atoms of ground (resp. upper) 
states |1), and |3) (resp. |2)) with the respective energies 
u>x, 0J3 and u>2- (see Fig. I). The pulses are delayed at 
the entrance of the medium such that the pump pulse is 
switched on first. This pulse sequence is referred to as 
intuitive sequence with respect to the well-known coun- 
terintuitive order in stimulated Raman adiabatic passage 
(STIRAP) processes We assume pulse durations to 
be much shorter than the relaxation times, such that 
losses from the upper state, and the decoherence between 
the two ground states due to collisions and laser phase 
fluctuations are negligible. The effect of the dissipation is 
considered in Section IV. C. For an exact two-photon res- 
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FIG. 1: Schematic diagram of the three-level atomic system. 



onance, the corresponding Hamiltonian in the resonant 
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approximation (RWA) reads in the basis {|1), |2), |3)} 





H= I -0, 



-o; o 

A -fi, 



-fi* 



(1) 







with fi PjS = n PyS E PjS /2H = |J7p )S |e l¥,p ' s the Rabi frequen- 
cies of the pump and Stokes laser fields, ^ PjS the corre 
sponding dipole moments, A = lj 2 — u)% — oj p the one- 
photon detuning. The instantaneous eigenstates write 
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with O 2 = | rip | 2 + \£l s \ 2 the generalized Rabi frequency, 
and the mixing angles defined as tan# = |f2 p /f2 s |, < 
9 < tt/2 and tan 20 = 20/A, < ip < tt/2 (0 < ip < 
7r/4 if A > 0). Note that, although for one atom the 
Rabi frequencies can be taken as real, their phases change 
in general during the propagation and become complex. 
The corresponding eigenvalues read 



Aft 2 = 

A d = 



i (A - + A 2 



0. 



= —Cltanip, (3a) 

= £1 cotan-0, (3b) 
(3c) 



State |1) is connected to the bright state \bi) (I&2)) 
when ip — v (ip — ► tt/2), corresponding to A > (A < 
0), and 6 — > tt/2 which is satisfied at early times when 
the pump pulse is switched on first. We will consider 
without loss of generality A > 0. The bright state is 
up to a phase the state solution during the interaction if 
the adiabatic conditions arc fulfilled: 



|A&! — Ad I » \[i0 — (ifp — (p s ) cos0 sin#] cosi/'l 



(4a) 



Aft, — A/, 



» 



iip + -(y> p sin 2 8 + ip s cos 2 9) sin 2-0 , (4b) 



where the dot corresponds to the derivative with respect 
to time. This is satisfied for 



|AT| > 1, fl 2 T/\A\ ~ |AT|^ 2 > 1 



(5) 



with T the time of interaction and ip, f2 considered during 
the pulse overlapping. The first inequality means that the 
spectral width of the pulses should be much smaller than 
the one-photon detuning, and the second one that this 
width should be much smaller than the Stark shift of the 
levels. We have obtained the latter inequality considering 
the case of interest fl < | A|, which leads to sin-0 ~ ip and 
cos ip ~ 1. Note that the population of the upper state is 
proportional to ip 2 and in order to reduce the losses from 
this level (of rate T) one should require additionally 



TTip < 1. 



(6) 



The Maxwell equations in the running coordinates 

T) = X, T = t — x/c (7) 

and in the slowly varying amplitudes approximation read 



dVt p 
drj 



iq p a x a 2 , 



dr] 



iq s a 3 a 2 



(8) 



where q p , s — 2iru> PtS fj 1 2 s N/hc are the coupling coeffi- 
cients, with N the atom density in the medium, and ai 
the atomic population amplitudes determined from the 
Schrodinger equation 



d 



H 



(9) 



with <p = [ox a 2 03]* (t denotes the transpose). For 
simplicity, we consider in what follows equal oscillator 
strengths q p — q s = q. 



III. EFFECTIVE PROPAGATION EQUATIONS 

A. First order equations 

Combining the Schrodinger equation §§§ and the prop- 
agation equations J5J) leads to 



dn 2 

drj 



9 I I 

OT 



drj 
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da^a\ 
dr 



(10) 



Using the adiabatic solution of the Schrodinger equa- 
tion, we get finally the following effective propagation 
equations for the angles 9, ip and the relative phase 

ip = ip p ~ ip s : 
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(11a) 
(lib) 
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This system of equations is valid only when the deriva- 
tives at all orders of ip, 9 and <p exist and are small such 
that the adiabatic approximation is satisfied. This ex- 
cludes in particular field envelopes of bounded domains 
for which there exist discontinuities of the derivative at 
a certain order. 

We remark that this system of equations (jTTJ) contains 
the time derivative of 9 and ip, and can be thus inter- 
preted as taking into account the first non-adiabatic cor- 
rections. We indeed recover these equations starting from 
the Maxwell equations pH|) to which we insert the solu- 
tion to the Schrodinger equation including the first order 
non adiabatic corrections, i.e. keeping linear terms in d T 9 
and d T ip. This calculation can be also interpreted as an 
adiabatic evolution, not in the usual adiabatic basis, but 
in the first order superadiabatic basis. 
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B. Solutions 

Analytical solutions to (fTTj) c an be found by the stan- 
dard characteristic method fill ]: 

1>{ti, t) = MO, 0(r), t) = 9 (0, <p(ri, t) = <po(t), (12) 

where ^o(C) = <P(v = 0, r = C), 0o(O = = 0,r = 
and <po(0 = V 3 ( 7 ? = 0; T — £) are the boundary conditions 
given at the entrance of the medium, and C = C( ? 7? r ) and 
£ = £(77, t) are solutions of the respective equations 



C = ^-^cos 3 2^o(C), (13a) 
ng(i)di = wcos 4 Vo(C)(2 - cos 2^o (0) (13b) 



C 



with f2(ry, r) = 0(77 = 0, r = t) = £lo(t) the generalized 
Rabi frequency given at the entrance of the medium. The 
pump and Stokes Rabi frequencies are respectively deter- 
mined from Q p — rising and fl s — ilcos(9. 

Through the definition tan 2ip = 20/ A, Eq. (jllaj) can 
be interpreted as describing the dynamics of the general- 
ized Rabi frequency (for a constant A) Cl(r], r) = fio(C)- 
We conclude that it propagates in the medium through 
the non-linear time £ and with the non-linear velocity v 
such that 1/v = l/c + gcos 3 2^o/A 2 , smaller than the 
light velocity c, and the delay in the medium for small 
angles ip at a certain length L is equal to 



qL 



qL 



T m = L(l/v - l/c) = —r cos A 2ip ~ 72 



(14) 



From Eq. fl3b]) . the non-linear time £ is always larger 
than £, meaning that the mixing angle propagates in the 
medium with the velocity greater than the light velocity 
[l3L [l4l ] . This can be also inferred from Eq. (jllbj) , where 
the non linear velocity appears negative. 

The solution (|12p shows that, if at the medium en- 
trance the relative phase ip of the pulses is constant, it 
remains constant during the propagation in the medium. 



C. Limitations for the adiabatic passage 

We show several limitations of solution (Tj"2")) . all re- 
lated to the adiabaticity of the process. They can be 
interpreted in terms of energy leading to the maximum 
propagation length [see Eq. (|16p ] and to the maximum 
time [see Eq. (TITl) ]. Non-adiabatic phenomena that lead 
to a singularity in the solution, and corresponding to a 
reshaping of the pulses, are prevented through the limita- 
tion pTj) [or ([22]) in the limit of small angle ip] . Reshaping 
of the pulses is also prevented by ([25]) . 



1. Energy 

At the beginning of the interaction r — > —00, i.e. Q — > 
-co according to Eq. (|13a[) . the quantity £(rj, r — > -co) 



is estimated by the equation 

c£(t),t-> -00) 



n (r)dT ~ gry. 



(15) 



We infer that £(t],t —* —00) is a monotonic increasing 
function of r\ changing from — 00 at r\ = to +00 at 
T] — Vmux, where rj max is estimated from 



+00 



nl(r)dT ~ qrj n 



(16) 



This gives a relation between the energy of the pulses and 
the maximum theoretical propagation length. The physi- 
cal interpretation of the condition (fTB)) for the length can 
be expressed in terms of energy: the number of atoms 
whose population can be fully transferred by adiabatic 
passage from state 1 1) to state |3) in the medium cannot 
exceed the number of photons in the pump and Stokes 
fields. 



2. Time 



The definition of £ (|13a|) shows an additional limita- 
tion: For a given length 77, one can get a large value for 
C(?7,t) when taking a large r, i.e. when t 3> qr)/A 2 . 
The definition of £ (|13bj) shows then that for a too large 
C, £ does not exist. For a given 77, the maximum value 
£ — > +00 is obtained for T max estimated from 



nl(t)dt ~ qr\. 



(17) 



Tmax has to be interpreted as a maximum time until 
which the adiabaticity is preserved. Beyond this value 
of T max , one expects non-adiabatic effects (beyond the 
first superadiabatic order). Eq. (TTT]) means that, during 
the propagation corresponding to larger 77, the maximum 
time from which adiabaticity is broken becomes smaller. 



3. Reshaping of the pulses 

Adiabatic conditions ([5]) have to be revised as follows 
when considering propagation. The time derivatives in 
Eq. (in the running coordinates) become: 

dtp _ dtpodC_ 86_ _ d0od^ dp _ dfo_d^ 

dr ~ dC, cV 8t ~ d£ 0r' cV ~ d£ St' [ ' 

This shows that the adiabatic conditions are the ones at 
the entrance of the medium provided that 



9c <i, SUi. 



(19) 



0T 0T 

The breaking of these conditions corresponds generally 
to the formation of shock-waves during the propagation. 
We obtain for the variable C 

dC _ 1 

07 ~~ 1 - (6qri/A 2 )(dip /d() cos 2 2^o sin 2^' 



(20) 
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At large propagation lengths the denominator in the 
right-hand side of this expression can become very small, 
that would lead to a strong increase of the derivative of 
the left-hand side. Thus, in order that the adiabaticity 
be preserved during propagation it is necessary to add to 
conditions j5]): 



6"0o tan 2-0Q 



T 



a I t. o / 1 L COs3 ^ i 
60 o tan2^ o -^ ~ < 1- 



(21) 



We infer that for relatively small angles ip and propaga- 
tion lengths L at which the group delay of the general- 
ized Rabi frequency is of the order of the characteristic 
interaction time, the adiabaticity in the medium does not 
break down. Considering for simplicity the case of small 
angles ipo, we get the condition at the length r\ = L: 



qTL < 



(AT) 2 
12 



0„ 



(22) 



with O max the maximum value of O. At a given peak 
Rabi frequency, a larger A allows thus one to increase 
the medium length for which pulses can propagate with- 
out reshaping. One can interpret this effect with the 
argument that a larger detuning weakens the effective 
interaction of the pulses with the medium. 
With respect to the variable £, one obtains 



dj _ o|(o oc 

dr fig(f) dr 



6qrj dtp 2 
~~ A 5 " ~d£ C0S ^° Sm ^ ( 



(23) 



Taking into account the preceding condition (|2"Tj) , this 
leads to 



3£ ng(Q d{ 

dr ~ dr' 



(24) 



This shows that the part 8 (and also <p) of the nonadia- 
batic coupling at the entrance of the medium is multiplied 
by the scale factor Qq(£)/Qq(£) during the propagation. 
The adiabaticity will break down when this ratio is too 
large, i.e. when £ > £ for late times. In practice, one can 
estimate the condition on the propagation length rj = L 
by preventing this situation, imposing £ — £ -C T, which 
gives from Eq. (|13b[) 



qTL < (O max T) 5 



(25) 



Inequalities (|2"Tj) and (f2"5|) guarantee the adiabaticity con- 
ditions during the pulse propagation, preventing a signif- 
icant reshaping of the pulses. 



IV. POPULATION TRANSFER BY B-STIRAP 

We consider the process of complete population trans- 
fer from state 1 to state 3 by adiabatic passage through 
the bright state using a so-called intuitive sequence of 
pulses, i.e. with the pump field switched on before 
the Stokes field with a delay > 0. The aim of 



this section is devoted to the theoretical interpretation 
of this process. We consider pulses of Gaussian en- 



velopes: O p 



Or, 



; exp[-((i + r d /2)/T p )% O s 



O s , max exph((<-T ti /2)/T s ) 2 j 

The main limitation for an adiabatic transfer for a sin- 
gle atom is given by ([5]). It gives (i) lower and upper 
limits for the detuning (at given fields) and (ii) a better 
adiabaticity for a larger pulse overlapping (that can be 
achieved with larger Rabi frequencies for given delay and 
pulse shapes). As studied above, there are four additional 
limitations given by [fig]), (JT7J), (J5J) [or |J52J) for small an- 
gle ipo] and (|25|) . The energetic argument (jTHJ) confirms 
the need of large Rabi frequencies to preserve the adia- 
baticity during the propagation. The condition (fT7|) on 
the other hand favours faster processes, i.e. finishing be- 
fore r max . We show below that this is the case for larger 
detunings. Condition (I22|l confirms that a larger detun- 
ing allows the propagation over a longer distance. Condi- 
tion ([23)1 is shown below to give the practical limitation 
on the medium length for which nearly complete popu- 
lation transfer can occur. Due to the interaction with 
the atoms, the pulses change their shapes. This corre- 
sponds to the violation of the adiabaticity conditions, as 
is studied below. We consider numerical calculations in 
different situations and interpret them in terms of the 
limitations described above. 



A. Application in experimental conditions 



We first apply the analysis to the experiment presented 
in [||, for which P;max T s = 108.6, Q s , max T s = 110.5, 
r d /T s = 1.4, T p /T s = 0.8, and AT S = 50. Figure dis- 
plays the dynamics of population transfer to state 3 for 
various lengths qT s x of the medium and the correspond- 
ing projections of the state vector \<fi) on the eigenstates. 
The analytic solution (fT2|) fits well the numerics until 
the complete transfer fails. The efficiency of the pop- 
ulation transfer decreases as the pulses propagate into 
the medium. The population transfer occurs completely 
in the medium up to qT s x w 4. For qT s x = 5, one 
can already notice a partial transfer of population. To 
achieve an effective population transfer one should pro- 
vide an adiabatic interaction to make state vector \<p) fol- 
low \bi) state and a smooth change of the mixing angle 
from 7r/2 to 0. However, during the propagation, 
one notices that state vector \<p) acquires a component 
along the dark state |d) leading to some final population 
to state |1). This is clearly seen from the right column 
of Fig. 2 which shows the instantaneous projections of 
state vector \4>) onto the dressed states. We can see that 
at qT s x = 2, there is no projections onto states I&2) and 
while at qT s x — 5 a noticeable projection onto state 
\d) appears. As a result the mixing angle changes 
and no more satisfies the conditions that are necessary to 
realize the transfer process. Indeed, in Fig. 3 we present 
the time evolution of the mixing angle for the same 
parameters as in Fig. 2. One can see that already at 
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Final state population State vector projections 




-2 2 -3 -2 -1 1 2 3 




-2 2 -3 -2 -1 1 2 3 




-3 -2 -1 1 2 3 -2 2 



Time (units of T s ) 



FIG. 2: (Color online) For (a) qT s x = 0, (b) qT 3 x = 2, and 
(c) qT s x — 5, as a function of the running time t — t — 
x/c, left column: Population transfer to state 3 determined 
numerically (full line) and from the analytical solution (|12|l 
(dotted line); right column: Projections |(d|</>)| 2 (dashed line) 
and | (6 2 1 0) | 2 (dotted line). 




-2 2 

Time (units of T s ) 



FIG. 3: Time evolution of the mixing angle 9 for qT s x = 
(dashed line) and qT s x — 5 (solid line). The inset shows an 
enlargement of the values of the mixing angle near the end of 
the interaction. 



propagation length qT s x = 5 the final value of is 
different from 0. 




-3-2-10123 

Time (units of T) 

FIG. 4: (Color online) Dynamics of the population transfer 
Pj = IO !0)| 2 (lower frame) corresponding to the pulses (upper 
frame) at the entrance of the medium x = for TA = 1000 
and TO, 

,max — 

TO 

s,max — 100. 



B. Interpretation and improvement of the transfer 

To interpret the maximal length leading to the com- 
plete transfer and to show how to improve it, we analyze 
the transfer taking similar parameters of the experiment. 
We take for simplicity equal peak Rabi frequencies and 
equal durations for both pulses. 

The main limitation is here given by Eq. (|17[) . It 
shows indeed that a process corresponding to a popula- 
tion transfer lasting too long, i.e. beyond r max , leads to 
an inefficient non-adiabatic transfer when propagation is 
considered. To illustrate this effect, we consider below 
two limiting cases with the same Rabi frequencies and 
pulse durations (fi p , max = f2 s , max , T p = T s — T), the 
delay Ta = 1.3T, and satisfying (J5J): one with a large de- 
tuning defined as A > !] max (see Fig. [4J, leading to a 
negligible population in the intermediate state, and an- 
other one with a small detuning defined as A ~ f2 max (see 
Fig. [5]), leading to a noticeable population in the inter- 
mediate state. At the entrance of the medium, r\ — 0, the 
population transfer for the large detuning case is accom- 
plished approximately during the overlap of the pulses 
(see Fig. 0]), while, for the small detuning, it occurs for 
the duration of both pulses (see Fig. [5]). (The calcu- 
lation with the small detuning corresponds to conditions 
that are very similar to the experiment described above.) 
Thus condition (|17[) is easier to satisfy for a larger de- 
tuning. The effect that the transient population in the 
intermediate state is smaller for a larger A corresponds 
to a weaker interaction of the pulses with the medium, 
and the condition (f2"Tj) [or (|22p] is also easier to satisfy. 
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One thus expects an adiabatic transfer over much longer 
propagation length for the case of large detuning. This 
is confirmed by the numerics. Figure [6] shows that, al- 
ready for qTx = 7, the population transfer is not com- 
plete. One can see that the corresponding analytic solu- 
tion fits well the numerics (until it stops at r « 2.35T), 
before the Stokes pulse is off and the transfer is complete. 
This time corresponds with a quite good approximation 
to the one determined from the estimation (fT7|) giving 
T max ~ 2.33T. Figure [B] shows that the nonadiabatic cor- 
rections correspond to a projection onto the dark state 
that grows roughly monotonically as a function of time, 
reaching its peak value near r max . 



100 




-1 1 

Time (units of T) 

FIG. 5: (Color online) Same as Fig. Hbut for TA = 50 and 

Fr2p ;max — FT2 s , ma x — 100. 

For the large detuning case, the limitation of the pop- 
ulation transfer during the propagation is given by (|25[) . 
Figure [7| shows the dynamics for qTx m&x = 100 which 
corresponds roughly to the maximum length with an effi- 
cient population transfer (it reaches here approximately 
99%). With respect to the condition (|25j) . this corre- 
sponds to gTL/(f2 max T) 2 = 0.01. Inspection of the ob- 
tained pulse shows that this non-adiabatic effect occurs 
near the end of the pump pulse, where it is significantly 
reshaped featuring a longer tail. Such a reshaping breaks 
the initial intuitive sequence of the pulses since when the 
Stokes falls down the pump pulse is not any more com- 
pletely switched off. Figure [7] shows that the nonadi- 
abatic corrections correspond to a projection onto the 
dark state. The maximal length obtained for the large 
detuning case is well beyond the one obtained for the 
small detuning case as discussed above. 

We conclude that in order to increase the propagation 
length for given identical pump and Stokes peak ampli- 
tudes, besides taking large detuning A [in the limit of 
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-1 1 

Time (units of T) 

FIG. 6: (Color online) Dynamics of the population trans- 
fer Pj = \(M)\ 2 (middle frame) for TA = 50, ™ p , max = 
Tfi s , max = 100 and at qTx — 7 (corresponding to Fig. [5] 
for x = 0); the absolute value of the corresponding Rabi fre- 
quencies (upper frame); and its corresponding non-adiabatic 
projections onto the dark state Pd = |(dj<^)| 2 and the bright 
state Pb 2 = |(&2|0)| • The dotted lines for the populations are 
determined from the analytic solution (|12l) . (|13[) . 



Eq. J5J], we should take the delay between the pulses as 
large as possible, while keeping an overlapping between 
the two Rabi frequencies of sufficiently large area (much 
bigger than one). 

We remark that the obtained maximum lengths of the 
medium where the complete population transfer occurs 
are here well below the lengths that would lead to a sin- 
gularity in the pulse shapes, as given by condition (|22p . 



C. Practical implementation and examples 

Equations J5} and ((6]) are the conditions that must be 
met in order to achieve an efficient population transfer 
by b-STIRAP in a medium. We consider below vari- 
ous media and experimental conditions that satisfy these 
conditions. 

In the gas phase, the limitations imposed by the losses 
from the three-state system are mainly due to the inho- 
mogeneous linewidth of the excited state [through Eq. 
([6])] and to the ionization by the fields, mainly from 
the excited state. In the case of a warm vapor, the 
inhomogeneous linewidth is mainly determined by the 
Doppler broadening, typically V — 500MHz for a warm 
alkali-metal-atomic vapor. In a regime corresponding to 

~ A, leading to a relatively short propagation length 
and a quite large transient population in the excited 
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Time (units of T) 

FIG. 7: (Color online) Same as Fig. [6]but for TA = 1000 and 
Tfip, m ax = Tfi s ,max = 100 and at qTx = 100 (corresponding 
to Fig. |4]for x = 0). The small discrepancy between numerics 
and the analytic solution is due to nonadiabatic corrections 
of order higher than one. 



state as typically shown in Figs. 6 and 9, all the cri- 
teria are fulfilled for a pulse duration of order T = 40 
ps. In the conditions of Fig. 9, we have fl — A/2, 
(£!/A) 2 rT = 0.005 (corresponding to a negligible loss), 
AT = 40, n 2 T/A = 10 and Q = 500 GHz. This im- 
poses to find a particular system with a limited ioniza- 
tion from the excited state with such a Rabi frequency, 
which corresponds to typical intensities of 30 MW/cm 2 
(for a transition strength of order of one Debye). For 
a long propagation length regime (as typically shown in 
Fig. 7, for f2/A = 0.1) with a small transient popula- 
tion in the excited state, one obtains T = 1 ns as the 
typical pulse duration. In the conditions of Fig. 7, this 
gives (fl/A) 2 TT = 0.005, AT = 1000, and Q 2 T/A = 100 
corresponding to the Rabi frequency fl = 100 GHz, i.e. 
to a generally weakly ionizing field of 1 MW/cm 2 (for a 
transition strength of order of one Debye). We remark 
that such a situation leads to TT ~ 0.5. However it is not 
expected to give a large loss due to the small transient 
population in the upper state (of order (Jl/A) 2 ~ 0.01). 
We have checked numerically this assumption as shown 
in Fig. 8, where various loss rates are studied. The loss 
rate TT — 0.1 does not show a noticeable loss in the pro- 
cess, while the loss rate FT = 0.5 gives a small loss of 
less than 3 percents. 

The situation is less restrictive in cold medium (such as 
Bose Einstein Condensates system) where T is bounded 
by the natural rate of losses from the excited state which 
is of order of 10 MHz. In this case the pulse duration 
can be increased to T ~ 2 ns for the short propagation 



Time (units of T) 



FIG. 8: (Color online) Population transfer in the conditions of 
Fig. [7] but for various loss rates: TT = (full line), TT = 0.1 
(dashed line, almost undistinguishable from the full line) and 
TT = 0.5 (dotted line). 



length regime case (A ~ O) and to T ~ 50 ns for the long 
propagation length regime CI/ A ~ 0.1. In these cases, the 
Rabi frequencies range from £1 ~ 2 GHz to il ~ 20 GHz 
which are well below the ionization regime in general. 

For solids, one of the main additional difficulties comes 
from their large optical inhomogeneous broadening. For 
example, in rare-earth-doped crystals, that are prefer- 
able within the context of coherent optical behavior for 
/ — / transitions, the upper state inhomogeneous broad- 
ening is of the order of 10 GHz. However, there are some 
techniques like hole-burning techniques that allow one 
to select a subset of the ions within a particular spectral 
range and thus to reduce the large inhomogeneous broad- 
ening. This allows one to increase the pulse durations up 
to the microsecond regime. This technique has been used 
in the experiment of Ref. Q that we have taken to apply 
our theoretical results (see Sec. IV A). There is usually 
another issue due to the level splitting As in "atom-like" 
systems (doped solids, color centers, quantum dots) in 
the range of 10 MHz - 10 GHz. This gives limits for the 
Rabi frequencies to address single levels: £1 <C As- 



V. PULSE DYNAMICS AND RESHAPING 

We here show the dynamics of the intuitive pulse se- 
quence propagating in the medium over long distance 
beyond adiabatic limitations. We choose parameters 
that show a reshaping already for small distances in the 
medium: tt max T = 20, r d /T = 1, and AT = 40. 

Figure [5] shows a typical dynamics: After a certain 
length there occurs a noticeable reshaping corresponding 
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Time (units of T) 



FIG. 9: (Color online) Time evolution of the normalized Rabi 
frequencies of the probe and Stokes pulses for different prop- 
agation lengths z = qTx: a) z = 0; b) z = 10; c) z = 20; d) z 
= 30; e) z=40; f) z=50. 

to a lengthening of the shape of the pump intensity at 
the tail leading to additional peaks (here one peak is 
shown). Already for z = 10, the population transfer is 
not anymore complete in this situation. The adiabaticity 
is strongly violated beyond this value. The leading edge 
of the stokes pulse is next also reshaped also showing one 
additional peak. The sequence of pulses is not any more 
intuitive at large distances. 



one-photon detuning weakens the adiabaticity of the in- 
teraction. 

The second difference concerns the derived criteria of 
adiabaticity in a medium. During the propagation of 
the two pulses there is an exchange of energy between 
them that causes their reshaping which breaks down the 
adiabaticity. For a counterintuitive pulse sequence, in 
the case of equal oscillation strengths, the adiabaticity is 
preserved during propagation at any propagation length. 
Only in the case of unequal oscillator strengths it can 
break down, mainly due to a singularity corresponding 
to a steepness occurring in the shape of the pulse (see ya, 
12]). In the case of an intuitive pulse sequence as shown 
in the present work the adiabaticity is broken down in a 
medium even for equal oscillator strengths. 

As noticed very recently [H, [3] , an intuitive sequence 
of the switching of the pulse features a reshaping that 
can lead to a super luminal propagation [l5|, [lg, [13, 03 • 
Our analytic solution should provide new insights about 
this effect and its limitation. Such a study is in progress. 

b-STIRAP can be viewed as an alternative popula- 
tion transfer method. The superluminal property of b- 
STIRAP can find applications. For instance, the simul- 
taneous use of the b-STIRAP and STIRAP processes in 
multi-level systems (e.g. tripod, M-system) can be an ef- 
ficient method for creating time controlled excitations of 
target states. We also anticipate that a combination of 
STIRAP/b-STIRAP is of interest in quantum informa- 
tion. The counterintuitive STIRAP scheme indeed allows 
the population to be transferred from 1 to 3 - while the 
same pulse sequence permits transfer back from 3 to 1 in 
an intuitive b-STIRAP scheme. This could find applica- 
tions for the implementation of quantum gates. 



VI. CONCLUSION AND DISCUSSION 

We have analyzed the complete population transfer by 
adiabatic passage when the pulses propagate in a medium 
in an intuitive sequence. Using the analytic solution, 
determined from the first order nonadiabatic corrections, 
we have derived limitations on the maximal length of the 
medium where the complete population transfer can take 
place. Assuming the adiabatic conditions fulfilled at the 
entrance of the medium [see conditions ([5])]. the complete 
population occurs over a longer distance in the medium 
(i) for a larger one-photon detuning (at given peak Rabi 
frequencies of the pulses) and (ii) for a larger peak Rabi 
frequency of the pulses (at a given ratio of the one-photon 
detuning with the peak Rabi frequency of the pulses) . In 
the limit A 3> fi, the limitation of the transfer appears 
as non-adiabatic corrections and a reshaping of the pump 
pulse through the condition ([23)) . 

We can list the differences between the techniques of 
STIRAP and b-STIRAP. Unlike STIRAP, b-STIRAP de- 
pends on the one-photon detuning. Increasing the one- 
photon detuning leads to an improvement of the prop- 
agation length at which the population transfer is com- 
plete in a medium. On the other hand, increasing the 
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APPENDIX A: METHOD OF 
CHARACTERISTICS FOR iji AND 6 

In this appendix we derive the solution of Eqs. (jTTJ) of 
the form 

|| + a«.)^=0, (Al.) 
|-*<=0, (Alb) 
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using the method of characteristics [ll|. We first solve 
Eq. (jAla.p using a characteristic parametrization r; = 
7](s) and t = t(s) such that 

diji 

ds 
drj 

ds 
ch 

ds 

ip is constant along the characteristics: ip(r),r) = 

-0(?/(O), t(0)). From (j A2b[) . we get r](s) = s choosing 
77(0) = 0. This entails 

fl>{tj,T)=i>{0,C) = MC), C = r(0) (A3) 
with ( = C(?7, t) solution of Eq. R2c|) : 

( = T-r 1 a(M0)- (A4) 



0. 


(A2a) 


1, 


(A2b) 




(A2c) 



This leads to Eq. (fT3a) . 

We next solve Eq. (|Alb[) rewritten as 



dd_ _ q(0 o )+6(^o) 89 
dz l + za>^ )ip'o(C)d( 



= 



(A5) 



using the change of variables 

z = r), ( = t - T]a(ip ) 



(A6) 



and denoting a'(ipo) = da/dipo, 0o(C) = dipo/dC- The 
method of characteristics allows one to solve Eq. (|A5[) : 



with z(s) = s and 

dC = a(V>o) +b(^ ) 
ds l + S a'(0oK(C)' 

The solution of the latter equation reads 



x exp 

which gives Eq. (I13b|) 



c a[MC)]+b[MC')} 

V>o(C) 



a'(^) 



(A7) 



(A8) 



(A9) 
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